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A  THERMAL  STRESS  SOLUTION  FOR 
MULTILAYERED  COMPOSITE  CYLINDERS 

Mark  0.  Withaiall 


ABSTRACT 

An  elrftic  thermal  stress  solution  is  presented  for  a 
multilayered  composite  c^inder  subjected  to  a  radial 
temperature  field.  The  solution  assumes  each  layer  to  be  at 
a  uniform  temperamre,  therefore,  many  layers  are  required 
to  model  regions  with  large  gradients  in  temperature. 
Previously  developed  closed*form  solutions  for  a 
monolayered  cylinder  subjected  to  a  uniform  temperature 
change  and  loadings  of  internal  and  external  pressures  are 
used  to  construct  the  multilayered  solution.  Results 
produced  by  the  theoretical  solution  have  been  compared  to 
those  generated  from  a  finite  element  analysis  with  excellent 
agreement.  Finally,  by  combining  a  previously  developed 
multilayered  stress  solution  for  loadings  of  internal  pressure, 
extern^  pressure,  and  axial  force  with  the  plane-strain 
thermal  stress  solution  of  this  report,  the  multilayered 
thermal  stress  solution  for  generalized  plane-strain  boundary 
conditions  is  shown  to  be  easily  obuined. 

NOMEMCLATURP 
a  -  inner  radius  of  a  layer 
Ay  •  elements  of  the  compliance  matrix  for 
an  orthotropic  material 
b  -  outer  radius  of  a  layer 

p  •  intemal  pressure  on  a  layer 

q  •  external  pressure  on  a  layer 

r  •  radial  position  within  a  layer 
a  •  therm^  expansion  coefficient 

AT  •  uniform  temperature  change  of  layer 

«  -  strain 

a  •  stress 

Subscripts 

a,b  •  evaluated  at  inner  and  outer  radial  location 

tj9^  •  radial,  hoop,  and  axial  directions 

Superscripts 

pq  •  identifies  stresses  and  strains  that  are  associated 
with  the  monolayered  stress  solution  for  intemal 
pressure  ’p’  and  external  pressure  ’q’. 

T  •  identifies  stresses  and  strains  that  are  as'ociated 
with  the  monolayered  stress  solution  for  thermal 
loading  AT. 


INTRODUCTION 

In  certain  cylindrical  pressure  vessel  applications,  such  as 
cannon,  there  exists  a  severe  thermal  environment  in  which 
the  vessel  is  exposed  to  high  temperatures  and  large  heat 
inputs.  In  recent  years  there  has  been  increased  interest  in 
exploiting  the  unique  qualities  of  composite  materials  in  the 
design  of  large  caliber  cannons.  Most  often,  compiosite 
materials  are  employed  as  jackets  in  a  multilayered  cannon 
construction.  In  order  for  design  and  analysis  of  these  types 
of  vessels  to  be  carried  out,  a  thermal  stress  solution  for 
thick-walled  multilayered  composite  cylinders  must  be  used. 
The  problem  of  constructing  the  multfiayered  thermal  stress 
solution  involves  applying  the  proper  boundary  conditions  to 
the  appropriate  monolayered  solutions.  The  two  monolayered 
solutions  used  iii  this  construction  process  involve  constant 
temperature  thermal  loading  and  loading  of  intemal  and 
external  pressures.  The  solution  involving  intemal  and 
external  pressure  is  necessary  to  preserve  continuity  of  the 
structure  as  each  layer  expands  or  contracts  due  to  its 
temperature  change.  In  bofo  of  these  stress  solutions,  each 
layer  is  considered  to  be  an  orthotropic  material  and  under 
plane-nrain  boundary  conditions,  i.e.,  each  layer  has  zero 
axial  strain.  The  monolayered  thermal  stress  solution  comes 
from  previous  work  by  Gerstle  and  Reuter  (1975)  of  Sandia 
Laboratories.  Their  solution  is  actually  a  special  case  of  that 
given  by  Tauchert  (1974).  An  interesting  and  important 
result  of  this  solution  is  that-unlike  isotropic  cylinders,  which 
are  stress-firee  fora  uniform  temperature  change -orthotropic 
cylinders  sustain  radial,  hoop,  and  axial  stresses  as  a  result  of 
a  uniform  temperature  change.  The  monolayered  stress 
solution  for  loadings  of  intemal  and  external  pressures  is 
given  I7  Lekhnitskii  (1963).  Both  the  work  of  Gerstle  and 
Reuter  and  that  of  Lekhnitskii  represent  closed-form 
solutions  for  the  loadings  and  boundary  conditions 
mentioned.  The  essential  components  of  these  solutions  are 
more  fully  discussed  below. 

CIE9METRY,  MATERIAL,  AND  LOADING  DEFINITION 
A  multilayered  cylinder  can  be  viewed  as  an  assembly  of 
many  single-layered  cylinders.  It  is  appropriate,  therefore,  to 
begin  with  a  review  of  the  monolayered  orthotropic  cylinder 
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I 


proUeni  for  thennal  loading  and  loadings  of  internal  and 
external  pressures.  It  should  be  mentioned  at  this  point  that 
the  nomenclature  used  in  this  report  has  been  chosen  to  be 
consistent  with  a  previous  work  by  Witherell  (1992) 
regarding  multilayered  stress  solutions. 

In  the  monolayered  case,  the  cylinder  is  assumed  to  be 
long  with  ends  that  are  fixed.  The  cylinder  has  an  inner 
radius  ’a’  and  an  outer  radius  V  (see  Figure  1).  The 


Fg.  1.  Monolayered  cylinder  geometry. 

cylinder  is  subjected  to  an  internal  pressure  ’p*.  an  external 
pressure  ’q’.  and  a  temperature  rise  of  AT.  /in  axial  force 
Fzi  is  necessary  to  enforce  the  fixed  end  constraint 
generated  by  the  pressures  ’p’  and  ’q’,  as  shown  in  Figure  2. 
Similarly,  an  axial  force  Fz4  is  needed  to  enforce  the  fixed 
end  constraint  for  the  temperature  rise  AT  (see  Figure  3), 
A  cylindrically-orthotropic  material  is  assumed  with  its 
principal  axes  coincident  with  the  cylindrical  coordinate 
qrstem  defining  the  cylinder  geometry.  In  constructing 
multilayered  composite  cylinders,  layers  are  most  commonly 
used  m  -f  and  -  fiber  wrap  angle  pairs,  where  each  pair  can 
be  viewed  as  a  sin^e  orthotropic  layer.  An  orthotropic 
material  is  characterized  by  twelve  independent  thermo* 
mechanical  material  constants  consisting  of  three 
engineering  moduli  (Ei.E2,E5),  three  shear  moduli  (G^,  Gg, 
G„),  three  Poisson’s  ratios  (vu,Va.V3,),  and  th.'ee  thermal 
expansion  coefficients  The  numben  IJZJ  indicate 

the  principal  material  directions  which  for  the  above 
assumptions  correspond  to  the  radial,  hoop,  and  axial 
directions  of  the  cylinder  In  addition,  since  the 

principal  material  directions  correspond  to  the  principal 
directions  of  both  the  cylinder  geometry  and  the  applied 
loadings,  shear  effects  are  eliminated,  and  the  number  of 
material  constants  necessary  for  the  analysis  reduces  to  nine. 


Fig.  2.  Internal  and  external  pressure  loading  of 
monolayered  cylinder  under  plane-strain 
boundary  conditions. 


Traction  Free  Inner 
and  Outer  Surfaces 


Fg.  3.  Uniform  temperature  change  loading 
for  monolayered  cylinder  under  plane- 
strain  boundary  conditions. 


MONOLAYERED  THERMAL  STRESS  SOLUTION 
The  equations  describing  the  thermal  stresses  for  an 
nnhotropic  cylinder  with  the  fixed  end  constraint  and 
traction-free  inner  and  outer  surfaces  are  given  by  Gerstle 
and  Reuter  (1975).  These  equations,  in  a  slightly  different 
form,  are  given  below.  A  superscript  T  is  used  to  signify 
thermal  loading. 
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Likewise,  tLe  hoop  strain  at  these  surfaces  is  given  by 


which  siniplifies  to 


eJ.-4.Ar 
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Similarly, 


which  simplifies  to 


eJ..4»Ar 


o'*-  -is®;* 
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‘Similar  to  what  was  done  to  find  Fz4,  F^i  can  be  found  by 
integrating  the  axial  stress  over  the  end  area  to  produce 


Again,  it  will  be  necessary  to  have  sttesses  evaluated  at  the 
inner  and  outer  surfaces  for  purposes  of  constructing  the 
multilayered  solution.  The  radial  and  hoop  stress  at  these 
two  locations  are  given  below 


oJJ  •  TAPp*TAQ-q 


4,-pj,Z)(*aC.‘-*-*(l-a)Cf'-l]*;,  (24) 


PLANE-STRAIN  MONOUYERED  SOLITTION  FOR 
INTERNAL  AND  EXTERNAL  PRESSURES 

The  equations  describing  the  stress  distribution  for  an 
oithotropic  cylinder  with  fixed  ends  and  loadings  of  internal 
and  external  pressures  arc  given  by  Lekhnitsk’i  (1963).  The 
form  of  these  equations,  given  below,  is  that  suggested  by 
O’Hara  (1987).  The  superscript  pq  is  used  to  signify  the 
solution  for  pressure  loading 


(26) 
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Fg.  4.  MSOL4  loading. 


PLANE.STRAIN  MULTILAYERED  THERMAL  STRESS 
SOLUTION  -  MSOU 

As  was  mentioned  earlier,  the  multilayered  thermal  stress 
solution  requires  both  the  monolayered  thennal  stress 
solution  and  the  monolayered  stress  solution  for  loadings  of 
internal  and  external  pressures.  This  is  because  as  an 
individual  layer  shrinks  or  expands  due  to  its  temperature 
change,  it  wU  either  push  or  pull  on  adjacent  layers.  This 
pushing  or  pulling  gives  rise  to  interface  pressures  at  the 
boundaries  of  adjacent  layers.  Figure  4  shows  the  loading 
for  the  MSOL4  solution.  In  order  for  continuity  of  the 
multilayered  cylinder  to  be  maintained  for  a  given 
temperature  distribution,  the  hoop  strain  of  two  layers  at 
their  interface  must  be  equal.  The  equation  defining  the 
hoop  strain  equivalence  at  the  interface  of  layer  i  and  i+l 
is  given  below 

•Uo  *  <fl(o  -  tZiM) 

By  substituting  the  appropriate  stress  componenu  from  the 
two  monolayered  solutions  outlined  above,  we  arrive  at  the 
following  equation: 


.*»(0A3TiT-  A.(i*I)ATO+l)*  ««-l)  ,3^, 

*Oaq(0*Ouq(t*l)~0 

where 

-Clma(0TAQ(0 


Ofl-  -WiaW-  ♦  fia(0-TAJX0  * 

and 

Ou-  (39) 

For  each  two*layer  interface  there  is  one  equation  defining 
the  hoop  strain  equivalence  condition.  For  a  c^inder  with 
TT  orthotropk  layers,  there  are  N>1  equations  and  N'l 
unknowns.  Setting  up  these  equations  in  matrix  form  and 
noting  that  p(l)  and  a(N)  are  the  prescribed  internal  and 
external  pressures  of  the  overall  c^inder  results  in 
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where 

T(0-  Aj(0AIT0-  a/<*l)A2T(*l)  (‘‘D 


In  a  condensed  form,  the  matrix  equation  above  can  be 
written  as 


(42) 


It  should  be  mentioned  that  p(l)  and  q(N)  should  be  set  to 
zero  if  the  objective  is  to  find  the  multilayered  thennal 
stress  solution  for  traction-free  inner  and  outer  surfaces. 

Die  system  of  equations  given  above  can  be  easily  solved 
to  deteimine  the  interface  pressure  vector  [Q].  By  recalling 
that  q(i-l)  w  p(i),  we  can  now  calculate  the  stress 
distribution  in  layer  i  by  using  AT.  p(i),  and  q(i)  as  input 
into  the  two  monolayeied  stress  solutions  dismissed  earlier. 

Each  layer  of  the  cylinder  has  two  axial  forces,  Fz.  and 
Fz4.  that  make  up  the  total  force  required  to  enforce  the 
fixed  end  conditions.  By  summing  each  of  these  forces  for 
the  entire  cylinder,  we  can  determine  the  total  axial  force 
(Fzu)  from  the  pressure  loads  and  the  total  axial  force  (Fj^) 
from  the  thermal  loading,  as  given  below 

F’b,  •  jF’afO 


THE  MULDLAYERED  THERMAL  STRESS  SOLUTION 
FOR  GENERAUZED  PLANE-STRAIN  BOUNDARY 
CONDIDONS  .  MSOLT 

In  the  multilayered  plane-strain  thermal  solution 
(MSOL4)  discussed  in  the  previous  section,  the  ends  of  the 
cylinder  are  fixed  and  the  axial  strain  is  zero  for  each  layer. 
Generalized  plane-strain  boundary  conditions  address  the 
case  in  srtiich  each  layer  of  the  cylinder  has  equivalent  but, 
in  general,  nonzero  axial  strain.  Both  the  fixed  ends  and  free 
ends  type  of  boundary  conditions  would  be  special  cases  of 
the  generalized  plane-strain  class  of  problems.  Axial  loading 
can  also  be  modeled  when  this  type  of  boundary  condition  is 
made  part  of  the  solutions  formulation.  ^  combining 
MSOL4with  anonthermal  generalized  plane-strain  solution 
(MSOLPQF)  by  Witherell  (1992),  the  multilayered  thermal 
stress  solution  for  generalized  plane-strain  boundary 
conditions  (MSOLT)  can  be  easily  constructed. 

The  generalized  plane-strain  multilayered  stress  solution 
(MSOLPQF)  is  for  loadings  of  intern^  pressure,  external 
pressure,  and  axial  force.  This  solution  is  construmed  from 
three  separate  multilayered  solutions,  MSOLl,  MSOL2,  and 
MSOL3.  The  MSOLT  solution  is  very  simply  obtained  from 
the  MSOLPQF  solution  by  substituting  MSOL4  for  MSOLl. 
Making  this  substitution  also  requires  adding  one  more  term 
(Fz«)  to  the  FZN  constant  used  in  the  MSOLPQF  solution. 
The  new  equation  for  F2^  is  given  as 


RESULTS 

The,  multilayered  thermal  stress  solution  MSOLT  was 
implemented  in  a  FORTRAN  program  and  verified  with  a 
finite  element  solution.  In  comparing  theoretical  predictions 
of  maximum  absolute  values  of  radial,  hoop,  and  axial  stress 
with  those  produced  using  finite  elements,  a  difference  of  less 
than  0.4  percent  was  found.  This  close  agreement  is  to  be 
expected,  since  the  finite  element  solution  will  approach  the 
MSOLT  solution  as  the  element  size  approaches  zero. 

An  example  of  a  solution  for  a  cylinder  with  an  internal 
radius  of  1  inch  (2J4  cm),  external  radius  of  2  inches  (S.08 
cm),  with  50  orthotropic  layers  is  shown  in  Figures  5, 6,  and 
7.  The  cylinder  b  under  generalized  plane-strain  boundary 
conditions  with  a  net  axial  force  of  zero.  The  layup  consists 
of  alternating  hoop  and  axial  layers  of  a  graphite/epoxy 
composite.  An  IM6  graphite  fiber  was  used,  and  the 
composite  had  a  60  percent  fiber  volume  ratio.  The  laminate 
material  properties  were  generated  using  an  average  of  four 
micromechanics  models  and  are  shown  in  Table  1. 

The  cylinder  is  subjected  to  a  linear  mdial  temperature 
distribution  starting  at  50*F  at  the  inner  layer,  decreasing  to 
1*F  at  the  outer  layer.  Also,  remember  that  the  temperature 
within  a  layer  is  uniform.  Figures  5  and  6  show  the  thermal 
stress  distribution  resulting  from  this  temperature  gradient. 
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TABLE  1.  UYER  MATERIAL  PROPERTIES  FOR  IM6/EPOXY 
69  PERCENT  FIBER-VOLUME  RATIO 


Fiber 

Direction 

(Mpsi) 

(lOr^rF)  1 

Er 

E, 

«, 

a. 

«. 

Hoop 

1.274 

25.40 

1.274 

0.0157 

0J127 

0.3956 

0.2304 

-0.0263 

0.2304 

Axial 

1.274 

1.274 

25.40 

0J956 

0.0157 

0J127 

0.2304 

0.2304 

■0.0263 

1  Mpsi  a  6.89  GPa 
IT*  a  9/5*C‘ 


RADIAL  POSITION  (In.)  RADIAL  POSITION  (in.) 


Fig.  5.  Axial  and  radial  thermal  stress 
distribution  (1  Ksi  a  6.89  MPa, 
t  in.  a  2J4  cm). 


Rg.  6.  Hoop  thermal  stress  distribution 

(1  ibi  a  6.89  MPa,  1  in.  a  2.54  cm). 


Because  the  inner  and  outer  surfaces  are  traction-ftee,  the 
radial  stress  at  these  locations  is  zero.  The  large 
fluctuations  in  hoop  and  axial  stresses  occur  at  axial-hoop 
layer  toundaries.  The  thermal  strain  distribution  is  shown 
in  Figure  7.-  The  radial  strain  curve  has  some  jaggedness  to 
it  because  of  jumps  in  temperature  from  layer  to  layer. 


Fig.  7.  Axial,  radial,  and  hoop  thermal 

strain  distribution  (1  in.  «  2.54  cm). 


CONCLb.’MONS 

Qosed-foim  elastic  stress  solutions  for  monclaycted 
cylinders  under  thermal  and  pressure  loading  were  used  to 
construct  the  plane-strain  multilayered  thermal  stress 
solution.  In  addition,  it  was  shown  that  this  solution  can  be 
simply  extended  to  the  case  of  generalized  plane-strain 
boundary  conditions.  The  solution  has  been  compared  with 
unite  element  results  with  excellent  agreement,  ^ally,  the 
multilayered  thermal  stress  solution  has  been  successfully 
integrated  with  a  finite  difference  heat  transfer  solution  to 
predict  transient  stress  distributions  in  multilayered  cannon. 
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